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Abstract 

We compute by QCD sum rules the matrix elements of the relevant four-quark 
operators appearing in the expression of the A b inclusive decay rates at the 
order l/m|. The results suggest that l/m| corrections are not responsible of 
the observed difference between the lifetime of and 
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1 Introduction 



An interesting problem of the present-day heavy quark physics is represented by the 
measured difference between the A& baryon and Bd meson lifetimes: t(A°) = 1.18±0.07 ps 
and t(B°) = 1.56 ± 0.05 ps [Q. As a matter of fact, the deviation from unity, at the 
level of 20%, of the ratio r(A b )/r( J B ci ): r(A°)/r( J B°) = 0.75 ± 0.05, is in contradiction 
with the naive expectation that, at the scale of the b quark mass, the spectator model 
should describe rather accurately the decays of the hadrons containing one heavy quark. 
The new measurement by the Delphi Collaboration at LEP of the average b— baryon 
lifetime: r{b — baryon) = 1.25 ± 0.11 ± 0.05 ps 0, together with the CDF result T(Afe) = 
1.33 ± 0.16 ± 0.07 ps [[|, although pointing towards a larger value of r(Ab)/r(_Bd), still 
confirms that this ratio is considerably different from unity. 

In principle, the ratio t{Ki ) )/t{B ( i) can be computed in QCD. As a matter of fact, 
a field theoretical approach has been developed for the analysis of the inclusive decay 
rates of the hadrons Hq containing one heavy quark Q 0, |!|. The method is based on 
an expansion in the inverse heavy quark mass mq, in the framework of the operator 
product expansion (OPE), and it has provided us with several hints on the hierarchy of 
the lifetimes of these hadronic systems. The calculation, however, involves the matrix 
elements of a number of high dimensional quark and gluon operators. For the D = 5 
operators defined below such matrix elements have been theoretically calculated, or can 
be inferred from the experimental measurements. As for the D = 6 operators, in the case 
of heavy mesons their matrix elements can be obtained invoking the factorization ansatz, 
and therefore it is possible to express them in terms of quantities such as the leptonic 
constant fs, whose estimates can be found in the literature. On the other hand, the heavy 
baryon matrix elements of the D = 6 operators cannot be obtained by factorization, and 
indeed they have been estimated only in the framework of constituent quark models, with 
uncertainties whose size is difficult to assess. Since the D = 6 A^ matrix elements may be 
responsible of the large difference between r(A b ) and r{B d ), it is interesting to compute 
them by field theoretical approaches; this paper is devoted to a calculation based on the 
method of QCD sum rules. 

Before reporting on this calculation, let us briefly summarize the main aspects of the 
QCD analysis of the inclusive decay widths of the heavy hadrons. The starting point is 
the transition operator T{Q — > Xf — > Q) j|: 

f = tjd 4 x T[C w (x)C ] w {<d)} (1) 
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describing an amplitude with the heavy quark Q having the same momentum in the initial 
and final state. Cw is the effective weak Lagrangian governing the decay Q — > Xf. The 
inclusive width the hadron Hq can be obtained by averaging T over Hq and taking the 
imaginary part of the forward matrix element: 

r( _ , = g Iru < W Q > 

^ m Hq 

The main idea to calculate the r.h.s of eq.(§) is to set up an operator product expansion 
for the transition operator T in terms of local operators Oi. 

f = J2 C i°i (3) 

i 

with Oi ordered according to their dimension, and the coefficients Cj containing appropri- 
ate inverse powers of the heavy quark mass rriQ. The lowest dimension operator appearing 
in (^|) is O3 = QQ. The next gauge and Lorentz invariant operator is the D = 5 chromo- 
magnetic operator C G : Oq = Q^cr^G^Q, whose hadronic matrix element 

i/c{Hq) = < H Q m^ QlHQ > (4) 

measures the coupling of the heavy quark spin to the spin of the light degrees of freedom 
in the hadron Hq, and therefore is responsible of the mass splitting between hadrons 
belonging to the same si multiplet (s£ is the total angular momentum of the light degrees of 
freedom in Hq). In the case of 6-flavoured hadrons this mass difference has been measured, 
both for mesons (M B * - M B = 42.0 ± 0.6 MeV 0) and S 6 baryons (M E » - M Sb = 
56 ± 16 MeV [§). 

The matrix element of QQ over Hq can be obtained using the heavy quark equation 
of motion, expanded in the heavy quark mass: 

QQ = g 7 °g + ^ - ^ + 0(m Q 3 ) ; (5) 

Ok is the kinetic energy operator O n = Q(iD) 2 Q whose matrix element 

, < H Q \Q(iD)*Q\H Q > 
^(H Q ) _ (6) 

measures the average squared momentum of the heavy quark inside Hq. On the other 
hand, the Hq matrix element of QrfQ is unity (modulo the covariant normalization of 
the states). 



The number of independent operators appearing in (|3]) increases if the 1/mg term is 
considered. Such operators can be identified in the four-quark operators of the type 



Ol = QTq qTQ (7) 

where T is an appropriate combination of Dirac and color matrices. 

In this way, a complete classification of the various contributions to the inclusive decay 
rates can be obtained for the different hadrons Hq. In the expression for T(Hq — > Xf): 

T(H Q - X f ) = Tl\Al + 4 + 4. + 1 ( 8 ) 

m Q m Q 

the A[ factors, that together with Tq depend on the final state Xf, include perturbative 
short-distance coefficients and nonperturbative hadronic matrix elements incorporating 
the long range dynamics. The partonic prediction for the width in (H) corresponds to the 



m; 
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leading term TP art {H Q -> X/) = T f A f , with A f = 1 + c f a s /7c + 0{a 2 s ) and T f 
differences among the widths of the hadrons Hq emerge at the next to leading order in 
1/rriQ and are related to the different value of the matrix elements of the operators (9, of 
dimension larger than three. 

It is important to notice the absence of the first order term ttiq 1 in @, a result obtained 
by Chay, Georgi and Grinstein ||, and Bigi, Uraltsev and Vainshtein |J. 

The occurrence of operators of the type in eq.(^) is an appealing feature of the expan- 
sion (||), as far as the determination of the inclusive widths is concerned. As a matter of 
fact, contrarily to the D = 5 operators Oq and O n which are spectator blind, the D = 6 
operators give different contributions when averaged over hadrons belonging to the same 
577 (3) light flavour multiplet, and therefore they are responsible of the different lifetime 
of, e.g., B~ and B s , Af, and H&. The spectator flavour dependence is related to the mech- 
anisms of weak scattering and Pauli interference 0, both suppressed by the factor m^ 3 
with respect to the parton decay rate. 

As for the differences in the lifetime of mesons and baryons, they could already arise 
at the order m^ 2 , due both to the chromomagnetic contribution and to the kinetic energy 
term in @. In particular, the kinetic energy term is responsible of the difference for 
systems where the chromomagnetic contribution vanishes, namely in the case of A& and 
Sft having the light degrees of freedom in S— wave. However, the results of a calculation 
of f/%. for mesons 0] and baryons |TIJ support the conjecture, put forward in |12| , that the 



kinetic energy operator has the same matrix element when computed on such hadronic 
systems. The approximate equality of the kinetic energy operator on Bd and A& can also 
be inferred by considering that, to the leading order in I/thq, /^(A b ) can be related to 
fi^(Bd) and to the heavy quark masses by the expression (which assumes the charm mass 
m c heavy enough for a meaningful expansion in l/m c ) JTc 



/4(A 6 ) - /4 (B d ) ~ - mbmc . [(M B + 3Mb* - 4M A J - (M D + 3M D * - 4M A J] . (9) 
2(m b - m c ) 

Using present data and the CDF measurement M\ b = 5623 ± 5 ± 4 MeV [Q (in |] the 
value 5639 ± 15 MeV is reported) eq. <g) gives /4(A 6 ) - f4{B d ) ~ 0.002 ± 0.024 GeV 2 , 
where the error mainly comes from the error on M^ b . The QCD sum rule outcome for 
hI(Hq) is fil(Bd) — /-4(Ab) ~ 0.6 GeV 2 , with an estimated uncertainty of about 30%. 
This result implies that the differences between meson and baryon lifetimes should occur 
at the itiq 3 level, thus involving the four-quark operators in eq.(|7]). They can be classified 



as follows [151: 



0\_ A = Ql1h<Il QlIuQl 

°s-p = QrIl QlQr 

T^_ A = Qhl^—qL QlI^—Ql 

n- P = qryk&y Qr (10) 

with q^i = ^~2^Q an d A a the Gell-Mann matrices. 

For mesons, the vacuum saturation approximation can be used to compute the matrix 
elements of the operators in ([3~0|) : 

<r R \<1« IR ^ ( m b + m q ^2 f'B q M k nn 

< B q \U v _ A \B q > VSA - \~^f B ) < B q\ U S-p\ B q>VSA- ^ (HJ 

< B q \Ty_ A \B q >vsa = < B q \Tg_ p \B q > V sa= . (12) 

Therefore, the matrix elements are expressed in terms of quantities such as Jb and the 
quark masses, and the resulting numerical values can be used in the calculation of the 
lifetimes, with the only caveat concerning the accuracy of the factorization approximation 
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The vacuum saturation approach cannot be employed for baryons; in this case a direct 
calculation of the matrix elements is required, for example using constituent quark models. 
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A simplification can be obtained for A&, as noticed in [TE\, using color and Fierz 
identities and introducing the operators 

6 q v- A = Qh.Qi ti-fd (is) 

and 

&s-p = QWr qIQr (14) 

(i and j are color indices). As a matter of fact, the A b matrix elements of the operators in 
fllOD can be expressed in terms of < A b \O v _ A \A b > and < Ab\Oy_ A \Ab >, modulo I/tuq 
corrections contributing to subleading terms in the expression for the inclusive widths. 

The matrix element of O v _ A and O v _ A can be parametrized as 

{ & v _ <^\Ov-a\M> _f 2 B M B . , 

and 

<A b \O q v _ A \A b > = -B <A b \O v _ A \A b > (16) 

with B = 1 in the valence quark approximation. 

For f B = 200 MeV and r = 1 eq. (|H>D corresponds to the value: (Oy_ A )A b = 4.4 x 
10~ 3 GeV 3 . In ref.iT^] the A c matrix element of O v _ A has been computed using a bag 
model and a nonrelativistic quark model; the results (O v _ a )a c — 0.75 x 10~ 3 GeV 3 and 
(O v _ a )a c — 2.5 x 10~ 3 GeV 3 , correspond to r ~ 0.2 and r ~ 0.6, respectively. An 
analysis using the model in [17j can also be found in ||. 



Larger values of the matrix elements have been advocated by Rosner using the 
values of the mass splitting — H b and S* — E c , and assuming that the A b and wave 
functions are similar: r ~ 0.9 ± 0.1 taking M|» — M| 6 = M|» — M| c , or r ~ 1.8 ± 0.5 
using the Delphi measurement in 0. 

It is worth supplementing the information from constituent quark models by estimates 
based on field theoretical approaches, for example QCD sum rules. As a matter of fact, 
a large value of r, namely r ~ 4 — 5, would explain the difference between r(A b ) 



and t(B c i) |L5" |. As we shall see, the application of the QCD sum rule method to the 
calculation of the matrix element of an operator of high dimension presents a number of 
disadvantages; nevertheless, interesting and quite reliable information can be obtained for 

< A b \Oy_ A \A b >. 
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2 QCD sum rule calculation of (Oy_ A )\ h 



A quantitative estimate of the matrix element (O v _ A )\ b can be obtained by the method 
of QCD sum rules [19|| applied to a suitable correlator in the heavy quark effective theory 
(HQET). Let us consider the three-point correlation function: 

U cd (uj,uj') = (l+^) CD n{u,u') =i 2 J dxdy < 0\T[J c {x)6 v _ A {0)J D (y))\0 > e M^)^'(vy) 

(17) 

of the spin ~ local fields J(x) and J(y) {C and D are Dirac indices) and of the operator 
Oy_ A in eq.(^). The variable uj (uj') is related to the residual momentum of the incoming 
(outgoing) baryonic current: 

p» = mbV » + k» (18) 

with fc M = ojv^. 

If the baryonic currents J and J have non-vanishing projection on the A& state 

<0\Jc\Mv)>= IaMc (19) 

(ip v is a spinor for a A& of four-velocity v) with the parameter f\ h representing the cou- 
pling of the current J to the A& state, the matrix element {O v _ A )^ b can be obtained 
by saturating the correlator ( |TTD with baryonic states, and considering the double pole 
contribution in the variables u and w': 

U^ d (u,u') = (6 v _ A )J-f- -i- + . . . (20) 

at the value uj = uj' = Aa 6 . The mass parameter Aa 6 represents the binding energy of 
the light degrees of freedom in A^ in the static color field generated by the b— quark: 

M Ab =m b + A Ab (21) 

and must be derived within the same QCD sum rule theoretical framework. 

A suitable interpolating field for A&, in the infinite heavy quark mass limit, has been 
proposed by Shuryak [2(J . It is set up by a function of the quark fields: 

J c {x) = e^{q T \x)TTqi{x)){h k v ) c {x) (22) 

where T means transpose, i,j and k are color indices, and C is the Dirac index of the 
effective heavy quark field h v (x) related to the Dirac field Q by: 

h v (x) = e im ^ x ^-Q(x) . (23) 
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The matrix r is the A& light flavour matrix: 

n 1 \ 

(24) 



JL ( 1 
I -1 



corresponding to zero isospin. As implied by the spectroscopy of baryons containing one 
heavy quark in the limit m& — > oo, in the A& the light diquark is in a relative + spin- 
parity state; this feature can be described by the current (p2"]) by two possible choices of 
the Dirac matrix T: 

r« = c 75 , r( 2 ) = C757 (25) 

(for a At at rest) where C is the charge conjugation matrix. In principle, the currents 
obtained using the matrices in ([25|), and a linear combination 

T = C l5 (l + bjf) (26) 



can be used in (|I7). As discussed in | IT| , there are arguments in favour of the choice 



b = 1 in (|26|); we shall come to this point below. 

Also the coupling f\ b can be derived by a sum rule, considering the two-point corre- 
lator: 

H CD {u) = (1 + ^)cdH(uj) =ijdx< 0\T[J c {x)J D {0)]\0 > (27) 
saturated by a set of baryonic states, in correspondence to the pole at u = A\ b : 

H Md (u;) = & 1 +... . (28) 
2 (A Af , - UJ) 

Therefore, the calculation of f\ h can be carried out following the same QCD sum rule 
approach in HQET; this analysis can be found in the literature [|T]], fH[]. 

Let us consider the correlators (pj]) and (|27j). In the Euclidean region, for negative 
values of u,u' the correlation functions (p7|) and ( p7|) can be computed in QCD, in terms 
of a perturbative contribution and of vacuum condensates. The results can be written in 
a dispersive form: 

IT OP >, J) = [ dado' Pn( ^ g/) , (29) 

H OPE {u) = [ da^\ (30) 
J (a — uj) 
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where possible subtraction terms have been omitted. The spectral function of ([T7|) read 
in HQET as 

Pn(a, a') = p% ert \a, a') + pg= 3 \a, a') < qq > + Pr ? =4) (a, «0 < —G 2 > 

71 

+ pg =5 \a, ( r')<qgcrGq>+p ( ° =6 \ ( T,CT')(<qq>) 2 + ... ; (31) 

a similar expression can be given for p#(cr). 

At the lowest order in a s the diagrams contributing to pn(cr, a') are depicted in 
fig.l. The perturbative contribution to the spectral function, obtained computing by the 
Cutkosky rule the imaginary part of the diagram in fig. la, has the following expression: 

ATM - 3^(1 + *)W - - a 4 + + (* ~ ✓» ■ (32) 

A comment is in order. The operator Oy_ A in the correlator ([17]) could give rise to 
non-spectator contributions through diagrams where two of the quark fields appearing 
in it are contracted in a tadpole. Such contributions, as noticed in j0|, do not affect 
the differences in lifetime for the various heavy hadrons, therefore they can be omitted 
assuming a normal ordering in the four-quark operators. Such contributions vanish when 
the non-perturbative spectral functions are computed. 

The contributions proportional to the quark condensate and to the mixed quark-gluon 
condensate can be derived by computing a class of diagrams of the type in fig.lb-d, with 
the result: 

p!f= 3 W') = -J^{9(a-o f )a' 2 (a-a') 2 + e(a'-a)a 3 (a-a') 2 + (T 2 a f2 } (33) 



b 

+ 9(a -a)(a' 2 + 9a 2 - 1(W) + 3<x 2 + 3a' 2 + W} . (34) 



Notice that in the adopted Fock-Schwinger gauge no gluon can be emitted from heavy 
quark leg in the infinite heavy quark mass limit. 

Terms proportional to the four-quark condensate come from diagrams of the type in 
fig.le, assuming factorization of the matrix element < qqqq >: 

" iD=6)f - ~" - ^rirKW + ^ 2( V) + ~ • (35) 

yb7r z 
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We do not include the gluon condensate contribution Pjj~^ (cr, a') since it represents the 
low momentum component of gluon exchange diagrams omitted in the calculation of the 
perturbative term. 

The inclusion of additional contributions proceeds in a similar way: we shall discuss 
in the following the consequences of the neglect of such terms. 

A sum rule for (Oy_ A )ji b can be derived by equating, according to the QCD sum rule 
strategy, the hadronic and the OPE representations of the correlator (|ITD . Moreover, 
invoking a global duality ansatz, the contribution of the higher resonances and of the 
continuum in n in ( pD| ) can be modeled as the QCD contribution outside the region 
(duality region) < uj < uj c , < uj' < u c , with uj c an effective threshold. Finally, the 
application of a double Borel transform to both the II and Yl had representation of the 
correlator in the momenta uj, uj': 

Bm— = -1 e~°/ E \ B{E X )-^— = JL e-*' E \ (36) 
cr — uj hi A — uj hi 

(and similar for uj' with the Borel parameter E2) allows us to remove the subtraction terms 
appearing in (|29|) , that are polynomials in the variables uj or uj' ( the Borel transform of 
a polynomial vanishes). Moreover, the convergence of the OPE is factorially improved by 
the transform, and the contribution of the low-lying resonances in ]J had is enhanced for 
low values of the Borel variables. The symmetry of the spectral functions in a, a' suggests 
the choice E\ — E2 — IE where E is the Borel parameter appearing in the QCD sum 
rule analysis of the two-point function (p7[). The final expression for the matrix element 
reads: 

^(l + 6) 2 exp(--^)(^_ A ) A6 = y o j dada'eM-^f)Pn(v,*') • (37) 



We use the standard values of the condensates appearing in (pi]): < qq >= (—240 MeV) 3 
and < qgoGq >= itLq < qq >, with = 0.8 GeV 2 |19| . The threshold parameter uj c 



has been fixed in the QCD sum rule determination of /a 6 and Aa 6 [11J|: uj c = 1.1 — 1.3 GeV. 



The parameter b appearing in the baryonic current J in (E3) has also been fixed in [J I 



studying the matrix element of the kinetic energy operator. The choice 6 = 1 allowed 
to obtain f Ab = (2.9 ± 0.5) x 10" 2 GeV 3 and A Ab = 0.9 ± 0.1 GeV. 

Using this set of parameters, we derive from eq.(|37D the result depicted in fig.2. A 

stability window is observed, starting at a value of the Borel variable E ~ 0.2 GeV and 

2 In the value < qq >— (—230 MeV) z has been used. The numerical results for f\ b and for Aa 6 , 
within the quoted uncertainty, are not affected by this choice. 
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continuing towards large values of E; in this range the result for {O v _ A ) Ai is independent 
of the external parameter E. It is known, however, that large values of the Borel variable 
do not provide us with interesting information, since in this region the sum rule is sensitive 
to the continuum model. Our sum rule (|37| ) is particularly affected by this problem, due 
to the high dimension of the spectral density Therefore, we are forced to consider a 
narrow region of E, close to the values E ~ 0.2 — 0.3 GeV; this is the Borel region already 
considered in the QCD sum rule analysis of f\ b and /i^(A&) [pjj . The variation of the sum 



rule result with E and with the continuum threshold uj c provides us with an estimate of 
the accuracy of the numerical outcome. We find: 

(O q v _ A ) Ab ~ (0.4 - 1.20) x 10- 3 GeV 3 , (38) 

a result corresponding to the parameter r in the range: r ~ 0.1 — 0.3. 



It is worth observing that the various contributions to the spectral function in 
have different signs and similar sizes, and that cancellations occur among the various 
terms. This is a common feature of the sum rule analyses of heavy baryon systems 
[ pi pi] , [?2j | , since the loop terms are suppressed by phase-space factors and are comparable 
with the nonperturbative corrections. A different procedure can be followed to soften 
this problem, using a partial resummation of the nonperturbative corrections through 
nonlocal condensates |f23| . Another possibility to test of the numerical result (|38| ) consists 



in assuming a local quark-hadron duality prescription [24j, that amounts to calculate the 



matrix elements of (O v _ A )^ b and f^ b by free quark states produced and annihilated by 
the baryonic currents in (|17D and (p7|), and then averaging on a duality interval in u,u'. 
The resulting equation for (O v _ A )^ b is simply given by: 

\ U V-A/A b ~ 1 rUc J—jr— , v\ ' 



20tt 



Jo c dcxcx 5 exp(— §) 



where in the denominator the perturbative spectral function of the two-point correlator in 
eq.(^) appears. It is possible to check the numerical outcome for the binding energy A\ b 
by this method: the result is reported in fig.3a-b where it is shown that the same value for 
such a parameter is obtained from 2 and 3-point sum rules. As for (O v _ A ) Ab , the result 
is depicted in fig. 3c; it corresponds to the value: (O v _ A )^ b ~ (0.6 — 1.2) x 10~ 3 GeV 3 in 
agreement with the result 



Let us consider, now, the parameter B in eq . (|16D . As it emerges considering the 
diagrams in fig. 1, in our computational scheme only valence quark processes are taken 
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into account, and therefore a sum rule for the matrix element in flltjD would produce the 
result B — 1. The calculation of the contribution corresponding to the diagrams of the 
type in fig. 1 immediately confirms this conclusion. 



3 Conclusions 

Within the uncertainties of the method, we have obtained by QCD sum rules small values 
for the matrix elements of (O v _ A )\ b and (Oy_ A )^ b , comparable with the outcome of 
constituent quark models. Results smaller then the value obtained using the mass splitting 
J^l — Sj, suggest that an experimental confirmation of this mass difference is required. 

From our results we conclude that the inclusion of I/tuq terms in the expression of 
the inclusive widths does not solve the puzzle represented by the difference between t(A&) 



and r{Bd). As a matter of fact, using the formulae in [15] for the lifetime ratio, the value 
in eq. (|38|) together with B = 1 gives: 

T{A b )/r{B d ) > 0.94 . (40) 

It seems to us unlikely that the next order contribution m^ 4 can solve the problem. If 
the measurement of r(A&) and T(Bd) will be confirmed in future, we feel that a reanalysis of 



the problem will be required, as suggested for example in |25[ . Meanwhile, it is interesting 
that new data are now available for other b— flavoured hadrons, e.g. Hf, |J, although with 
errors too large to perform a meaningful comparison with A&. Such new information will 
be of paramount importance for the complete study of the problem of the beauty hadron 
lifetimes, an argument that definitely deserves further investigation. 
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FIGURE CAPTIONS 



Fig. 1 

Set of diagrams contributing to the spectral function (|3l|). Diagram contributing to the 
perturbative term (a), to the D = 3 term (b), to the D = 5 term (c,d), to the D = 6 term 
(e). The thick lines correspond to the 6-quark propagator in HQET. 

Fig. 2 

Sum rule fl3~7|) for the matrix element {O v _ A )\ b as a function of the Borel variable E. The 
curves refer to the threshold parameter u c = 1.1 GeV (continuous line), uj c = 1.2 GeV 
(dashed line), uj c = 1.3 GeV (dotted line). 

Fig. 3 

Sum rule ( p9[ ) for the matrix element (Oy_ A )\ b obtained using the local duality ansatz. 
In (a) and (b) the mass parameter A, obtained from the two-point sum rule and from 
the sum rule for the four quark operators is depicted. The curves refer to u c = 1.1 GeV 
(continuous line), uj c = 1.2 GeV (dashed line), u c = 1.3 GeV (dotted line). 
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Fig.l 
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Fig.l 
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